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\ Abstract. Let Lu := V • (o(i)Vu) = — 5(x — y) in R 3 , < c\ < a(x) < C2, a(x) is a 

ff) 1 piecewise-smooth function with the discontinuity surface S which is smooth. It is proved that 

in a neighborhood of S the behavior of the function u is given by the formula: 

J (^a+r^r-y 1 +bR~ 1 }, y 3 >0, 
' u(x, y) = < , (*) 

' 1 (477a-)- 1 [r^ - bR- 1 ], y 3 < 0. 

Here the local coordinate system is chosen in which the origin lies on S, the plane X3 = is 
CNj ' tangent to S, a+(a_) is the limiting value of a(x) on S from the half-space X3 > 0, (x 3 < 0), 

r xy ■= \x-y\, R:= ^ p 1 + (\x 3 \ + |j/ 3 |) 2 , P ■= \J {xi ~ Vl) 2 + (X2 ~ V2) 2 , b := (a+ -a_)/(a + + 
a_). If S is the plane x 3 = and a(x) = a_|_ in x 3 > 0, a(a;) = a_ in 3:3 < 0, then (*) is the 
global formula for a in R 3 . Inequality for the fundamental solution for small and large \x — y\ 
follows from formula (*). 



1. INTRODUCTION. 

There are many papers on the behavior, as x — > y, of the fundamental solutions to the 
elliptic equations of the form 

' n du 

C$ ■ Lu := 22 ®3 [ a ij( x ) u j( x ^y)] = -$( x ~ v) m u j '■= -k — = djU (1.1) 

• • 1 uXj 

1,3 = 1 ■> 

for smooth coefficients a^-. Methods of pseudo-differential operators theory give expansion 
in smoothness of the solution to (1). In [LSW] existence of the unique solution to (1) with 
the properties 

< Cl r 2 - n < u(x, y) < c 2 r 2 ~ n , u £ Hl oc (R n \y), r := \x - y\, (1.2) 

is obtained under the assumption that a^- are bounded real-valued measurable functions 
such that 

n n n 

ai^^tf — a ij( x )tih — a 2 ^2 ' a ii °2 = const > 0. (1-3) 

i=l i,j = l i — 1 

Our purpose is to give an analytical formula for the fundamental solution of the basic 
model operator (1.1), namely the operator with 

aij(x) = 5 ij a(x), d(x) = { + ' 3 > ' (1.4) 

[a-, x 3 <0. 
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Here a + and a_ are positive constants, 

and u(x, y) is the unique solution of the problem: 

n 

Lu := d i (a(x)ui) = -5(x - y) in R n , (1.5) 

i=l 

[tt]|s = 0, [a(s)«jv] |s = 0, (1.6) 
the symbol [u]|s denotes the jump of u across S, that is, 

[u] = u + — U- ,u± := lim u(s ± sN) , 

£^0 

s is a point on S, N is the unit normal to S directed along x 3 , un is the normal derivative 
on S, S is the plane x 3 = 0, [ciun] '■= a+ufj — a-U~^. 

Problem (1.5)-(1.6) is important in many applications and is called a transmission prob- 
lem. The solution to (1.5)-(1.6) is sought in the class Hl oc (R n \y) n W^{R n ), where 
H 1 := W 1 ' 2 and W{* is the Sobolev space of functions whose distributional derivatives 
up to the order i belong to Lf oc . If 

f af, x 3 > 0, 
a ij {x) = { « (1.4') 
[ dij, x 3 < 0, 

and the constant matrices af- are positive definite, then there exists an orthogonal coordinate 
transformation which reduces af- to 5^ and a~- to XjSij, \j > 0. We do not give the formula 
for u(x, y) in this more general case. 

Finally note that for discontinuous coefficients equation (1.5) is understood in the weak 
sense, namely as the identity: 

/ a(x)ui(x,y)4>i(x)dx = (f>(y), (1.7) 

The identity (6) for u G Ht oc (R n \y) n W^(R n ) implies conditions (1.6). 

Let n = 3. The formula for the solution to problem (1.4)-(1.6), or the equivalent problem 
(1.4), (1.7) is given in Theorem 1.1. 

Theorem 1.1. The unique solution to problem (1.4)-(1.6) is: 

, . (^[H^ Vs >0, b:= a ^-, 
u(x,y) = < (1.8) 
I [?-£]> 2/3 < 0, r:=\x-y\, 

where R := V ( x i = Vi) 2 + i x 2 ~ 2/2) 2 + (\x 3 \ + I2/3I) 2 • 
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Corollary 1.1. The following inequality holds for r — > 0: 

\u(x,y)\ < c\x - y| _1 , (1.9) 

where the constant c > does not depend on x and y. 

Thus, the fundamental solution of the equation (1.1) with discontinuous senior coefficients 
has a different representation than the fundamental equation for the similar operator with 
continuous coefficients, but satisfies similar inequality for small \x — y\. 

A formula, similar to (1.8) can be derived by the same method for n > 3 as well. Formula 
(1.8) allows one to get asymptotics of u(x,y) and of V x u(x,y) as \x — y\ — > 0. Such 
asymptotics are useful in the study of inverse problems for discontinuous media [3]. 

In section 2 we prove Theorem 1.1. In section 3 various generalizations and applications 
are discussed. 



2. PROOF OF THEOREM 1. 



The proof is given for n = 3, but it holds with obvious small changes for n > 3. 

The idea of the proof is to take the Fourier transform of equation (1.5) with respect to the 
variables x := (x\,X2), to solve the resulting problem for an ordinary differential equation 
analytically, and then to Fourier-invert the solution of this problem. 

Let us go through the steps. 

Step 1. Let 2/1=2/2 = without loss of generality (since u is translation- invariant in the 
plane (xi,x 2 )). Denote 

w(Z,x 3 ,y 3 ) := e^ £ u{x,x 3 ;y)dx; u= — i— / we~^ x d£. (2.1) 

£ := (6,6), e = I?! 2 =e + & Denote w' := Jg. 

Let us Fourier-transform equation (1.5), with a(x) given in (1.4), and get 

w"(Z,X3,y 3 )-ew(Z,x 3 ,y 3 ) = { a +, (2.2) 

{ - — 6 ( X 3 ~ 2/3), X 3 < 0, 

w(6 +0,2/3) -«>(£, -0, y 3 ) = 0, a + w'(Z,+0,y 3 )-a-w'(Z,-0,y 3 ) = 0. (2.3) 

In what follows we omit £ in the variables of to, and write w(x 3 ,y 3 ) for brevity. Thus, to 
solves problem (2.2)-(2.3) and satisfies the condition 



w(±oo,y 3 ) = 
3 



(2.4) 



Assume that y 3 ^ 0. Then problem (2.2)-(2.4) has a solution and this solution is unique. 
A lengthy but straightforward calculation yields the formula for w: 

{ exp(-\£\x 3 -y 3 \) h exp[-\£\(\x 3 \ + \y 3 \)] n 

— mr + — + b mr + > f3 >o 
exp(-|g|^ 3 -y3l) _ h eM-\$\{\x3\ + \V3\)] . n 1 j 

2|£|a_ U 2|£|a_ ' #3 U 

where 6 := — — , a_,a + > 0. (2.6) 

a+ + a_ + v ' 

Step 2. The function u(x, y) is obtained from w by the second formula (2.1). Let us denote 
|£| := v, p := |ar| = \Jx\ + x 2 ,) and remember that 2/1=2/2 = 0. Since u; depends on |£| and 
does not depend on the angular variable, |£| := u, we have 



-1 /-OO /"Z7T -1 /-OO 

n = — — / di/i/ / e -™ pcosip wd(p = — dsvwJ Q (vp) (2.7) 



where Jo (a;) is the Bessel function and we have used the known formula: 

e ivpcosv d(p = J (vp). 
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We need another well-known formula: 



POO 



e-"* J {vp)dv = t>0 (2.8) 

VP + * 

From (2.5), (2.7) and (2.8) we get (1.8) with 2/1=2/2 = 0. Therefore, recalling the translation 
invariance of u in the horizontal directions, we get (1.8). 

Theorem 1.1 is proved. □ 
Remark 2.1 Note that the limits of u(x,y) as j/3 — > ±0 exist and are equal: 

n(x,^, + 0)=,(x^,-0) = 27rr(a+ 1 + Q _ ) . (2.9) 

A result similar to (2.9) is mentioned in [K, p. 318], however the argument [K] is not clear: 
the differentiation is done in the classical sense but the functions involved have no classical 
derivatives: they have a jump. 



3. GENERALIZATIONS, APPLICATIONS. 



This section contains some remarks. 
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Remark 3.1 First, note that if a(x) is a piecewise-smooth function with a smooth discon- 
tinuity surface S, s G S, a± = lim £ ^ a ( s =t e -^0> where N is the exterior normal to S at 
the point s, then the main term of the asymptotics of the fundamental solution u(x, y) in a 
neighborhood U s of the point s G S is given by formula (1.8) in which x, y G U s . This follows 
from the fact that the main term in smoothness of the solution to an elliptic equation in U s 
is the same as to the equation with constant coefficients which are limits of a(x) as x — ► s. 
In our case, this "frozen-coefficients" model problem is given by equations (1.4)-(1.6). This 
argument shows that the same conclusion holds if the coefficient a(x) in R 3 ^ and in R 3 is 
not smooth but just Lipschitz-continuous. 

Remark 3.2 In principle, our method for calculation of u(x,y) for the model problem 
(1.4)-(1.6) is applicable for the model problem (1.4') with anisotropic matrix. 
Remark 3.3 We only mention that our result concerning asymptotics of u(x, y) as \x — y\ — > 
for piecewise-smooth coefficients is applicable to inverse problems of geophysics and inverse 
scattering problems for acoustic and electromagnetic scattering by layered bodies. 
For example, if the governing equation is [R, p. 14]: 



k = const > 0, say k = 1, q(x) = 1 + p(x), where p(x) is a compactly supported real- valued 
function, p{x) G L^ oc (M 3 ), supp(p(x)) C R 3 . := {x : xs < 0}, a(x) = 1 +A(x), where A(x) is 
compactly supported piecewise-smooth function with finitely many closed compact smooth 
surfaces Sj C R 3 of discontinuity. Across these surfaces the transmission conditions (1.6) 
hold, and at infinity u satisfies the radiation condition. Then u(x, y) is uniquely determined. 

An inverse problem is: given g(x) and u(x, y) for all x, y G S := {x : X3 = 0} and a fixed 
k = 1, can one uniquely determine a(x), in particular, the discontinuity surfaces Sj? 

To explain how Theorem 1.1 can be used in this inverse problem, note that if two systems 
of surfaces Sj 1 ^ and and two functions a\ and 02 produce the same surface data on S 
for all x,y G S, then an orthogonality relation [R, pp. 65, 86] holds: 



where v(x) = 01 — 02, u m (x,y), m = 1,2, are the fundamental solutions corresponding to 
the obstacle D m (i.e., to a m and S^ m) ), D' 12 := R 3 \L>i 2 , D 12 := D 1 U D 2 . 

Let us prove, e.g., that dD 1 = dD 2 , using (3.1). If there is a part of dDi which lies outside 
D 2 , and s is a point at this part, then, assuming (for simplicity only) that v(x) is piecewise- 
constant and using for Vui and Vit2 formulas, which follow from (1.8) as y = z — > s, we 
conclude that the left-hand side of (3.1) is an integral which contains a part, unbounded 
as y = z — > s G dD\\ cf \x — y\~ 4 dx, c = const 7^ 0. This contradicts to (3.1). Therefore 
there is no part of dD\ which lies outside D 2 . Likewise, there is not part of dD 2 which lies 
outside dDi. Thus, dD\ = dD 2 . Similarly one proves that = for all j, provided 
(3.1) holds. 

A detailed presentation of such an argument is given in the paper by C. Athanasiadis, 
A.G. Ramm and I. Stratis, Inverse acoustic scattering by layered obstacle (in preparation). 

5 



V • [o(s)V«] + k 2 q(x)u = —5(x — y) 






(3.1) 
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